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We review the normal and superconducting state properties of the unconventional triplet superconductor 
Sr2RuC>4 with an emphasis on the analysis of the magnetic susceptibility and the role played by strong 
electronic correlations. In particular, we show that the magnetic activity arises from the itinerant electrons 
in the Ru d-orbitals and a strong magnetic anisotropy occurs (x -1 < X zz ) due to spin-orbit coupling. 
The latter results mainly from different values of the g-factor for the transverse and longitudinal compo- 
nents of the spin susceptibility (i.e. the matrix elements differ). Most importantly, this anisotropy and the 
presence of incommensurate antiferromagnetic and ferromagnetic fluctuations have strong consequences 
for the symmetry of the superconducting order parameter. In particular, reviewing spin fluctuation-induced 
Cooper-pairing scenario in application to Sr2Ru04 we show how p-wave Cooper-pairing with line nodes 
between neighboring Ru02 -planes may occur. 

We also discuss the open issues in Si'2Ru04 like the influence of magnetic and non-magnetic impurities 
on the superconducting and normal state of Sr2Ru04. It is clear that the physics of triplet superconductivity 
in Sr2RuC>4 is still far from being understood completely and remains to be analyzed more in more detail. 
It is of interest to apply the theory also to superconductivity in heavy-fermion systems exhibiting spin 
fluctuations. 
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1 Introduction 

The phenomenon of superconductivity remains one of the most interesting problems of condensed matter 
physics. In particular, in recent years the material science development has revealed several interesting 
systems where high transition temperature superconductivity was found, in particular the families of high- 
Tc cuprates Q] with a maximum T c of about 155 K. Soon after their discovery it was realized that the 
essential physics of cuprates takes place in the Q1O2 -planes which is believed to be responsible for the 
high transition temperature. 

In this connection the discovery of superconductivity in Sr2Ru04 with T c = 1.5 K |2| is of particular 
interest for several reasons. First, the crystal structure of Sr2Ru04 is identical to that of the parent com- 
pound of the high-T c superconductor La2Cu04 (see Fig. [Qfor illustration). Both kinds of materials are 
highly two-dimensional and as Fig. [Qshows the structure is almost identical to that of the L^-xSrzCuC^ 
superconductors. Both materials are oxides with conduction occurring in partially filled d— bands that are 
strongly hybridized with the oxygen p-orbitals. Therefore, it was generally believed that a comparison of 
the normal and superconducting properties of the cuprates and Sr2Ru04 will give a deeper understand- 
ing of the nature of the high-T c in the cuprates. However, it has been found that the differences between 
Si"2Ru04 and the cuprates are larger than their general similarities might suggest. In Sr2Ru04 super- 
conductivity occurs only at low temperatures and the normal state is a well-defined Fermi-liquid. This 
contrasts strongly with the anomalous normal state of the cuprates. Furthermore, it was soon found that the 
superconductivity in Sr2Ru04 is very interesting on its own. In particular, there are clear indications that 
the superconducting state is unconventional. For example, the transition temperature is highly sensitive to 
impurities 1 3 1 and nuclear quadrupole resonance(NQR) experiments do not show a Hebel-Slichter peak in 
the spin-lattice relaxation at T c 0. Shortly after the discovery of Sr2Ru04 it was suggested that super- 
conductivity might arise from odd-parity (spin-triplet) Cooper-pairs with total spin 5=1 and a non-zero 
angular momentum which is reminiscent of the phases of superfluid 3 He Q. The basis for this suggestion 
was the presence of ferromagnetism in the related compounds as SrRuC»3 and thus the expectation of fer- 
romagnetic fluctuations in metallic Si"2RuC»4. To support this picture the model phase diagram shown in 
Fig. |2]has been suggested. Here, one plots the phase of the ferromagnetic and superconducting members 




Fig. 1 (a) Crystal structures of (a) layered cuprate singlet superconductor La2- a; Sr I Cu04 and (b) layered perovskite 
structure of the triplet superconductor Sr2Ru04. Both structures are identical if La (Ba) is replaced by Sr and the 
Cu02-plane is replaced by Ru02-plane. 
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Fig. 2 Schematic phase diagram of the ferromagnetic and superconducting systems (Rodllesen-Popper crystals) 
Sr„ + iRu„03„+i, taken from Ref. (5). The number of layers is the parameter which determines the transition be- 
tween the two phases: MO=magnetically ordered and SC=superconducting. This phase digram suggests ferromagnetic 
excitations in the normal state of Sr2Ru04. 

of the so-called Rodllesen-Popper series (Sr n +iRu„03„ + i) as a function of the numbers of RuC>2-layers 
per unit cell, n. The infinite layer (SrRuC»3) is a ferromagnet with Tcurie ~ 165 K. For n=3 one finds 
Tcurie ^ 148 K and for n=2 the substance orders magnetically at Tcurie ~ 102 K. This demonstrates 
the tendency that Tcurie is reduced with decreasing layer number n and suggests that even for n=l, when 
supeconductivity occurs, one expects significant ferromagnetic fluctuations which may play an important 
role for triplet superconductivity in Sr2Ru04. 

Meanwhile, a number of experiments indeed point towards spin-triplet Cooper-pairing. The most con- 
vincing evidence comes from the 17 NMR Knight shift data which shows that the spin susceptibility 
is not affected by the superconducting state for a magnetic field parallel to the Ru02-plane j4). In Fig. 
|5]we show the corresponding experimental results. In conventional superconductors the spin part of the 
Knight shift measured by NMR decreases rapidly below T c due to the formation of singlet Cooper-pairs. 
On the other hand, in a triplet superconductor with S=l the spin part of the Knight shift should not change 
below T c for some field orientations, since the polarization induced by the weak external magnetic field 
and probed by NMR does not change. This behavior was observed in Sr2RuC>4 by Ishida et al. [4| and 
provides strong evidence for triplet Cooper-pairing. 

However, recently it has became clear that the situation is not that simple. For example, the inten- 
sive studies by means of inelastic neutron scattering (INS) reveal the presence of strong two-dimensional 
antiferromagnetic spin fluctuations at Q= (27r/3, 27r/3) and no sign of ferromagnetic fluctuations in the 
normal state of Sr 2 Ru04 [7 1. This fact, at first glance, cannot be compatible with triplet Cooper-pairing and 
thus the role of the antiferromagnetic fluctuations in the formation of superconductivity has to be clarified. 
Moreover, further analysis of the NMR data indicated the magnetic response to be strongly anisotropic 1 8 1. 
In particular, the measured longitudinal component of the spin susceptibility \zz is much larger than the 

transverse one %h , while from an isotropic model one expects xh = %Xzz- This anisotropy increases 

with decreasing temperature and reaches a maximum close to T c indicating its influence on the super- 
conducting properties. Therefore, in view of the importance of spin fluctuations for the unconventional 
superconductivity in Sr2RuC»4 their behavior should be understood in more detail. 

Another interesting question concerns the symmetry of the superconducting order parameter. While 
the analogy to the case of superfluid 3 He leads to the suggestion of nodeless p-wave superconductivity 
in Sr2RuC<4, recent analysis of the specific heat data |9| indicates a more complicated and anisotropic 
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Fig. 3 Results for the uniform spin susceptibility, %(q = 0,oj = 0) in the superconducting state of Sr2Ru04 as 
measured by NMR Knight shift |4|. One clearly sees that in contrast to the singlet Cooper-pairing (illustrated by the 
dashed curve) where spin susceptibility decreases upon cooling, the Knight shift in S^RuCk is unchanged by lowering 
T below T c . The Knight shift doe not decrease below T c , since the polarization induced by the external magnetic field 
does not change in the superconducting state in case of spin-triplet Cooper-pairs. 



behavior of the superconducting order parameter in the first Brillouin Zone (BZ). Therefore, the physics of 
Sr2RuC>4 is far from being completely understood. 

Note, an extensive review on the experimental situation can be found in Ref. [ilOjj. In contrast to that 
review, in the present review we mainly analyze the role played by strong electronic correlations and spin 
fluctuations in SraRuO^ So far their role in determining the superconducting and normal state properties 
has not been addressed in detail. Probably, this was due to the success of the conventional band theory in 
explaining the various properties of Sr2Ru04. However, its success is mainly due to the correct description 
of the Fermi surface topology in this compound. The Fermi surface shape, however, is hardly affected by 
the electronic correlations, while on the other hand the bandwidth and the energy dispersion are strongly 
modified. We will show that this is indeed the case for the ruthenates. Furthermore, we review the spin- 
fluctuation mediated Cooper-pairing scenario in application to Sr2Ru04 and address the important role of 
spin-orbit coupling in this system. In particular, we discuss the symmetry of the relevant superconducting 
order parameter within this scenario and, finally, investigate the role of impurities in Sr2Ru04. 



2 Superconductivity and magnetism in Sr 2 Ru0 4 

2. 1 Electronic structure and Fermi surface 

In Si"2RuC»4 the formal valence of the ruthenium ion is Ru 4+ . This leaves four electrons remaining in the 
4d-shell. Furthermore, the Ru ion sits at the center of a RuC»6 -octahedron and the crystal field of the O 2- 
ions splits the five 4d-states into threefold (t^g) and twofold (e s ) subshells as illustrated in Fig. |4] The 
negative charge of the 2 ~ ions causes the t2 g subshell to lie lower in energy. Note, that these orbitals 
(xz, yz, and xy) have lobes that point between the 2 ~ ions lying along the x, y, and z-axes. Electrons of 
these orbitals form the Fermi surface (the so-called a, f3, and 7-bands). We assume that the most important 
interaction of the carriers is with spin excitations described by the spin susceptibility x(q, w) and resulting 
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Fig. 4 Electronic structure of Sr2Ru04. The Ru 4+ ion corresponds to a 4d 4 level. The splitting of the e g and tig 
subshells is due to the RuC>6 crystal field. The orbitals d^, d xz , and d yz cross the Fermi level and form the a-, j3-, 
and 7-band. The resulting band structure is quasi-two-dimensional. Magnetic activity results from the tig subshell. 



from Ru a-, (3-, and 7-states. Then, we have 3 bands formed out of the oxygen p — d-ir hybridized orbitals 
and with a filling factor of 4/3. Since p — d-ir bonding is weaker than the p — d-a realized in cuprates, the 
admixture of the oxygen orbitals is smaller in ruthenates than in cuprates. For example, the contribution 
of the oxygen p-states close to the Fermi level is only 16 percent in Sr2Ru04. At the same time the 
contribution of the Ru 4d-orbitals is about 84 percent (11]. LDA calculations 1 12 13 1 indeed confirm the 
existence of the three bands crossing the Fermi level. 

The band structure of Sr2Ru04 is quasi-two-dimensional and the electronic dispersion along the c-axis 
is very small 1131 . This is because the highly planar structure of Sr2Ru04 prevents substantial energy 




Fig. 5 Calculated two-dimensional Fermi surface of S^RuCU using Eq. (1). The Fermi surface consists of three 
sheets: two with electron-like topology and one with hole-like topology. Note, the quasi-one-dimensional character of 
the q(xz)- and /?(yz)-bands. 

dispersion along the z-direction, due to the large interplanar separation of the RuC>6 octahedra. At the 
same time, in the afe-plane neighboring RuOe-octahedra share O ions which in turn are 7r-bonded to the 
Ru ions. Thus, the xy-orbital will form a two-dimensional band, while the xz- and yz-bands have only a 



8 



I. Eremin et al.: Unconventional superconductivity and magnetism in Sf2Ru04 



restricted one-dimensional character. Then, to describe the LDA results the following tight-binding energy 
dispersions can be introduced 



The lattice constants a = b have been set to unity. This model is a good parameterization of LDA calcula- 
tions 1141 if one chooses for d xy -, d zx -, and G?j, z -orbitals the values for the parameter set (e u , t x , t y , t') as 
(0.5, 0.42, 0.44, 0.14), (0.23, 0.31, 0.045, 0.01), and (0.23, 0.045, 0.31, 0.01) eV. The experimental investi- 
gation of the Fermi surface made by means of de-Haas-van- Alphen effect completely confirmed the LDA 
results 1151 . The resulting Fermi surface is shown in Fig. [5] Here, the a-band is hole-like centered around 
(tt, 7r) point of the Brillouin Zone, while (3 and 7-bands are electron-like ones centered around T-point. 
Furthermore, as one sees the a and /3-bands are quasi-one-dimensional and their degeneracy is removed 
by the introduction of the interband (between a and /3-bands) hopping t±= 0.025 eV I l*3l. 

Despite of the large difference in the magnitudes of the resistivity in the RuC»2 -plane and perpendicular 
to the planes, p c j p a b > 500, the temperature dependences of p c and p a b for T < 25 K both follow 
the Fermi-liquid behavior oc T 2 |3|. Above T — 25 K small deviations from this law occur. At higher 
temperatures p a b becomes linear with a temperature dependence similar to the cuprates. 

2.2 Dynamical spin susceptibility and electronic correlations 

As already mentioned in the introduction, the spin dynamics of Sr2Ru04 consists of ferromagnetic and 
incommensurate antiferromagnetic spin fluctuations at wave vector Qi = (2ir/3,2ir/3). Its role in the 
formation of triplet superconductivity in Si"2Ru04 requires special consideration. 

Long-range magnetic order is absent in Sr2Ru04. On the other hand the uniform magnetic suscepti- 
bility is much larger than the Pauli susceptibility of the non-interacting electrons. Indeed the 'ab-initio' 
calculations of the Stoner factor IN(0) (where / is the exchange parameter at q = 0, N(0) is the density 
of states at the Fermi level) gives the following values IN(Q) = 0.82 021, IN(0) = 0.89 (13). Using a 
simple so-called Stoner criterium, this indicates that Sr2Ru04 is in the vicinity of the instability to ferro- 
magnetism. In SrRuC»3 the corresponding Stoner parameter is IN(Q) — 1.23 1 16 1. Therefore, this system 
is a ferromagnetic metal with magnetic moment 1.6/zg and Curie temperature T c = 150 K. The reason for 
the enhanced ferromagnetic exchange fluctuations in SrRuC»3 in comparison to Sr2Ru04 is the stronger 
p — d hybridization and thus a larger contribution of the oxygen p-orbitals in the density of states at the 
Fermi level. 

However, what is interesting in Sr2Ru04 is not only the presence of ferromagnetic fluctuations at q = 
(from the Stoner instability), but also the incommensurate antiferromagnetic fluctuations at the wave vector 
Qi = (2-7r/3, 2ir/3) as observed in inelastic neutron scattering experiments |7|. What is the origin of these 
fluctuations? In fact they originate from the nesting properties of the quasi-one-dimensional xz- and yz- 
bands. In order to see this we show in Fig. [6]the results for the Lindhard response functions calculated for 
the different bands using their tight-binding dispersions 



where i refers to the band index and /j. is the Fermi distribution function. As seen in Fig. [6] due to the 
pronounced nesting of the xz and yz-bands their susceptibilities display peaks at Qi = (27r/3, 27r/3), 
while the icy-band does not show any significant feature. The response of the .Ty-band is enhanced due to 
the presence of the van-Hove singularity close to the Fermi level. Then, it becomes clear that the features 
observed by INS relate mainly to the magnetic response of the xz- and yz-bands. Furthermore, as has 
been shown in Ref. 1 17 1 that the Stoner enhancement related to the appearance of the incommensurate 
antiferromagnetic fluctuations is even stronger than the ferromagnetic one and is of the order of 1 . This 
means that Sr2Ru04 is almost unstable with respect to the spin density wave formation (SDW). Note, 



£o — 2t K cos k x 



— 2t y cos ky + 4t' cos k x cos k 



y 



(1) 




(2) 
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Fig. 6 Calculated real part of the Lindhard spin susceptibility xo (l refers to the band indices) in the normal state of 
Si"2Ru04 along the route (0, 0) —* (n, 0) — > (n, n) — > (0, 0) in the first BZ for the three different orbitals (xz,yz, and 
xy) crossing the Fermi level. Due to the nesting of the xz and yz-bands their susceptibilities show an enhancement 
at the incommensurate antiferromagnetic wave vector Qi = (2tt/3, 2n/3). The response of the xy-band is more 
isotropic, but significantly larger than in the normal metal due to the vicinity of the Van-Hove singularity. 



that both ferromagnetic and antiferromagnetic fluctuations are seen experimentally in the normal state of 
Sr2RuC>4 by means of 17 O Knight shift 1 18 1 and by inelastic neutron scattering |7|, respectively. 

These results suggest that an effective Hamiltonian for describing the physics in Si"2Ru04 is a two- 
dimensional three-band Hubbard Hamiltonian 

H = ^ ^ ek Q ak iQ(T ak,acr + ^ U a ni a ^n ia i , (3) 

where ak lQ(T is the Fourier transform of the annihilation operator for the d a orbital electrons (a — xy,yz, zx) 
and U a is an effective on-site Coulomb repulsion. The hopping integrals t^ a denote the energy dispersions 
of the tight-binding bands, see Eq.Q. 

As a result the spin dynamics in Sr2Ru04 is determined by the competition and the delicate balance 
between ferromagnetic and incommensurate antiferromagnetic fluctuations. The natural question arises 
what is the influence of these magnetic fluctuations on the triplet superconductivity in Sr2Ru04? In the 
following subsection we will study the formation of triplet of superconductivity in Sr2Ru04 driven by the 
spin fluctuation exchange. 

In the remainder of this section we discuss another aspect of the ruthenates, namely strong electronic 
correlations. As a matter of fact the success of 'ab-initio' methods in describing magnetism of Si"2RuC»4 
appear to invalidate, at first sight, the discussion of the role of strong electronic correlations. On the other 
hand, the effects of strong electronic correlations are small if the bandwidth of the conduction band W 
is much larger than the on-site Coulomb repulsion, U. Given that the bandwidth is of the order of 1 eV 
there is no reason to think that W >> U in Sr2Ru04. To estimate the value of U one has to remember 
that for 4d-electrons U has to be smaller than for 3d-electrons, since the average radius of the 4<i-shell is 
approximately two times larger than for the 3d-shell. Then a simple estimates gives ~ 1 eV. This 
means that Sr2Ru04 is in the regime of intermediate electronic correlations when U ~ W. 

As we mentioned earlier, despite of the importance of the electronic correlations in Sr2Ru04 the relative 
success of the band structure theory can be easily explained. The correct description of the ferromagnetic 
and antiferromagnetic fluctuations comes from the correct Fermi surface topology which is not affected by 
the electronic correlations. However, the bandwidth and the energy dispersion are strongly affected by the 
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electronic correlations in Si^RuOzi. For example, it has been shown by Liebsch and Lichtenstein 1 14 1 that 
in order to explain the bandwidth as measured by angle-resolved photoemission spectroscopy (ARPES) 
the electronic correlations have to be taken into account. In particular, using dynamical-mean-field theory 
(DMFT) they have shown that, while the Fermi surface topology is the same independent on electronic 
correlations, the bandwidth is different, namely the correlations reduce the bandwidth of the all bands in 
Sr2RuC>4. We will further address this issue discussion the phase diagram of Ca2- x Sr 2 ;Ru04. 



2.3 Superconducting state of S^RuC^: role of spin fluctuations and symmetry of the supercon- 
ducting order parameter 

Even before experimental results have indicated triplet superconductivity in S^RuO,}, its possibility had 
been theoretically predicted by Sigrist and Rice [5 1. In particular, they considered the triplet superconduc- 
tivity on a square lattice as an electronic analogy of the superfluid A-phase in 3 He. 

For the spin-triplet Cooper-pairing the wave function can be written as a matrix in spin space 

* = ffi(fc)|TT>+5a(*)(IU>+UT>) + sa(*)IU> = (J^j iffy) , (4) 

where the eigenvalues of the S 2 operator with projection +1,0, -1 have the form 

m>=(; ;) , iu>+ut>=(; j) , iu>=(s ;) . ® 

Another possibility to express the Cooper-pair wave function is to use the basis of the symmetric matrices 

itJGy — (ia x <Jy 7 i<Ty<Jy,i<J Z <Jy). (6) 

Then the Cooper-pair wave function has the form 

* = «(d(k)-a)=(-*« + ^ dx % dy ) . (7) 

The components of the vector d can be expressed linearly via the amplitudes of g a (k): 

9i = ~d x + idy, ,g 2 = dg, g 3 = d x + id y (8) 

Due to the Pauli principle the orbital part of the wave function with total 5=1 has to be odd, i.e. 
the orbital quantum number I — 1, 3, .... Therefore, one can speak about p, f, ... Cooper-pairing in a 
continuum system. In particular, the amplitudes g a (k) and the vector d(k) have to be the odd functions of 
the momentum k. 

In the tetragonal crystal structure such as that of S^RuCi, the order parameter d(k) has to belong to the 
corresponding irreducible representation of the group. There are several representations and 
the proper choice can be made by comparison to experiment. For example, the spin state of the Cooper- 
pair wave function can be obtained from the measuring the Knight shift in NMR and polarized inelastic 
neutron scattering experiments below T c . The orbital state can be measured by muon spin relaxation 
(/iSR) indicating time-reversal symmetry breaking. The square symmetry of the vortex lattice can also be 
interpreted as evidence for p-wave pairing. 

As we have shown in the introduction, the NMR Knight shift experiments indicate spin-triplet Cooper- 
pairing in Si"2Ru04 1 4 1. The same conclusion has been drawn from the polarized neutron scattering studies 
|20|. The time-reversal symmetry breaking detected by means of the /iSR |21 1 suggests the superconduct- 
ing order parameter E u with S z = and l z = ±1 (so-called chiralp-wave state): 

d z = zA (sin k x ± ismk v ), z = \S S = > (9) 
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Fig. 7 Lowest order RPA diagrams for the effective pairing interaction V e f / for singlet pairing resulting from the 
exchange of longitudinal and transverse spin and charge fluctuations. The solid lines refer to the one-particle Green's 
function and the dashed lines denote an effective Coulomb interaction U. The first diagram leads to a renormalized 
chemical potential. Note that for singlet pairing only an even number of bubble diagrams occur due to Pauli's principle. 



The amplitude of this gap is isotropic, Ak = Aoyk% + k^, and does not have a node at the Fermi surface 
of Sr 2 Ru0 4 . 

So far, the above analysis is quite general. Next, we assume the superconductivity is driven by the 
exchange of spin fluctuations. The interesting question in the ruthenates is why does triplet Cooper- 
pairing take place in the presence of strong incommensurate antiferromagnetic spin fluctuations at = 
(27r/3, 2ir/3) originating from the xz- and yz-bands and only relatively weak ferromagnetic fluctuations 
arising from the sy-band. In general, to answer this question the BCS gap equation (see Eq. ( I16» must be 
solved assuming different possible symmetries of the order parameter and corresponding pairing interac- 
tion for singlet, triplet with S z = projection and triplet with S z = ±1 projection Cooper-pairing. For the 
analysis of the interaction between quasiparticles and spin fluctuations (ferromagnetic and incommensurate 
antiferromagnetic) we use the generalized Eliashberg like theory for the interaction between quasiparticles 
and spin fluctuations \ 22\. We extend this theory to become a three-band theory, as necessary in Sr2RuC>4 
(23JEI. 

In the normal state both the self-energy and thermal Green's functions become matrices of 3 x 3 form, 
i.e. Gi.j ym (k, u> n ) and Ei.j m (k, ui n ), where i, j, m refer to the band indexes of the xy, yz, and xz-orbitals. 
The corresponding Dyson equation is given by 

~* -£)(k,w n ), (10) 

where m (k,oj n ) is the matrix of the bare Green's function determined via the tight-binding energy 
dispersions for the xy, yz, and xz-bands. The self-energy is given by 

S j;jVm (k,w„) = Vjj.mi^ Ml) x Gi tjim (k -q,u) n - m) , (11) 

q,l 

where m (q, uJi) is an effective interaction between quasiparticles and spin fluctuations. As shown in 
Fig. 7 it consists of an infinite series of diagrams including charge and spin fluctuations. This is similar to 
the case of cuprates, however, some important differences are present in the random phase approximation. 

Most importantly we must now consider the diagrams shown in Fig. [8] with an odd number of bubbles 
that contribute to the triplet pairing. This is in contrast to singlet pairing in cuprates where an even number 
of bubbles occur (see Fig. Q. Furthermore, the transverse (+-) and the longitudinal (zz) parts of the spin 
susceptibility are different. The reason for this (as we show later) is the presence of spin-orbit coupling. 



G(k,w„) 



12 



I, Eremin et al.: Unconventional superconductivity and magnetism in Sf2Ru04 



-OH + y-o-o-o- 



Fig. 8 Additional RPA diagrams of the lowest order for triplet Cooper-pairing with odd number of bubbles that refer 
to longitudinal charge and spin fluctuations. 

We discuss its role in detail in the next section. Then the effective pairing interaction in the 3x3 form 
including transverse and longitudinal spin fluctuations and also charge fluctuations is given by 

vy m ( q , Ul ) = ^X„?(q, w,) + ^;+-(q, + \Vit m {^ "l) , d2) 

where 

sp,zz sp,~\ — 

V sp ' zz = U 2 — X ° svzz , V sp ' + - = U 2 — ^9 — (13) 

i-u x s o p > i-UxT 

describe coupling to spin density fluctuations and 

ch 



ch _ jjl Xo 

i + u x f 



V = U \ rr ch ( 14 ) 



to charge density fluctuations, Xo h < XcT'" 1 are me irreducible parts of the charge and spin susceptibilities 
respectively. Note that in the Bethe-Salpether equation shown diagrammatically in Fig. f° r singlet 
Cooper-pairing it is necessary to have an even number of bubbles and include the ladder diagrams. In the 
case of triplet Cooper-pairing the contribution of the ladder diagrams is zero and an odd number bubbles 
diagrams occur due to Pauli's principle shown in Fig. [8] Since the Feynman rules require a (—1) for each 
loop an extra minus sign enters the gap equation via V t e . 

The magnetism in the ruthenates resulting from the quasiparticles in the t2 g -orbital is itinerant and thus 
the magnetic response is created by the same electrons that form also the Cooper-pairs. Then, for example, 
the irreducible part of the charge susceptibility \Q h is defined in terms of the electronic Green's functions 

xS' l (q) = ^E G ( k + c i) G ( k )' ( 15 ) 

k 

where G(k) is the single-electron Green's function (we omit here the band indices for simplicity). The 
longitudinal and transverse components of the spin susceptibilities are also calculated in terms of the elec- 
tronic Green's functions. In contrast to the cuprates, in the ruthenates Xo P ' ZZ an d Xo P H are different due 
to the magnetic anisotropy |4|. 

In the ruthenates the general matrix form of the superconducting gap equation for A(k, u>) is similar 
to the case of cuprates. To determine the superconducting transition temperature T c one must solve the 
following set of linearized gap equations to obtain A M (T): 

ApA^k,^) = -§ E E v Slm( k - k '' w » - ^) G U'^^i) 

k' ,u)j V ,m' 

x G mm /(-k', -Uj) A Mii /, m /(k', ujj). (16) 
Here, \ mu is an eigenvalue and T c is determined from the condition A(T C ) = 1, 
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Note, the interband coupling will provide a single T c for all three bands. The pairing potential 
(taken for singlet or triplet pairing) controls which state gives the lowest energy and determines singlet or 
triplet Cooper-pairing. Three possibilities may occur: 

(a) singlet pairing, 

(b) triplet pairing with the total spin S 2 = of the Cooper-pair 

wave function, 

(c) triplet pairing with the total spin S 2 = ±1. 



Note, a priori one cannot judge which pairing state is realized in the ruthenates due to the presence of 
both antiferromagnetic and ferromagnetic fluctuations. Therefore, one has to solve the gap equations for 
all three possibilities. Thus, for singlet pairing we use 

V? ) = lv£ + V+--±V eh . (17) 
For triplet pairing with S 2 = ±1 we take 

v Si = ~\v% ~ \Voh (18) 

resulting from the set of diagrams with odd number of bubbles (see Fig. |8j. Finally, for triplet pairing with 

S z = we get 

vgH = ~Vg-V+--±V ch . (19) 

Therefore due to magnetic anisotropy we have separated the longitudinal and transverse parts of the spin 
fluctuations. Note, transverse and longitudinal spin fluctuations contribute differently to the singlet and 
triplet Cooper-pairing. In particular, in the spin-triplet case the longitudinal and transverse spin fluctua- 
tions do not act additively as in the case of singlet Cooper-pairing. As a result, a strong reduction of the 
superconducting transition temperature T c is expected. The effect of charge fluctuations is expected to be 
small in the superconductivity in S^RuO,}. 

The pairing state symmetry is determined at T c by the eigenvector of Eq. ( 1 6) corresponding to A M (T) = 
1. As well as in the case of superfluid 3 He many other possible states besides the A-phase can occur. 
Note, that recent experiments indicate the occurrence of the line nodes in the superconducting state of 
Si"2RuC»4. These are the power rather an exponential laws of the specific heat, C(T) oc T 2 1251 . spin- 
lattice relaxation rate 1/Ti oc T 3 |26|, the heat capacity k oc T 2 |27|, and ultrasound attenuation |28|. 
Therefore, other symmetries of the superconducting order parameter like p-wave symmetry with nodes or 
/-wave symmetry which are other than simple nodeless p-wave symmetry of the superconducting gap have 
to be considered. Note, the largest eigenvalue of Eq. dl6l corresponds to a minimum in the free energy and 
thus will yield the symmetry of the superconducting order parameter A/ in S^RuOj. One could expect 
that the formation of the node should correspond to the influence of the incommensurate antiferromagnetic 
fluctuations. Let us see whether the simple nodeless p-wave symmetry is possible if the IAF at the wave 
vector Q=(27r/3, 27t/3) are present. 

Using appropriate symmetry representations 1 5 1 we discuss the solutions of the gap equation assuming 



p, d, or /-wave symmetry of the order parameter in the RuC>2 -plane |29|: 

d p (k) = Aoz(smk x a + ismk y a), (20) 

Ad(k) = A (cosfc x a — cos k y a), (21) 

d/ 2 2 (k) = Aqz(cos k x ci — cos k v a) (sin k x a + i sin k v a). (22) 



Here the z-unit vector refers to the d z component of the Cooper-pairs as observed in the experiment 1 32 1. 
We will present a possible explanation for this below. These symmetries of the superconducting order 



14 



I, Eremin et al.: Unconventional superconductivity and magnetism in Sr2Ru04 




Fig. 9 Symmetry analysis of the order parameter for triplet pairing in the first BZ for k z =0 (solid curves) and k z = 
■k /2c (dotted curves), a, f3, and 7 denote the FS of the corresponding hybridized bands. The wave vectors Qi and 
are the main wave vectors contributing to the susceptibility (without spin-orbit coupling). For f x 2_ y 2-wave symmetry 
the nodes of the real part of the order parameter are shown (dashed lines) and the regions + ( - ) where the / a .2_ H 2-wave 
superconducting gap is positive (negative). Note, for the real part of a p-wave order parameter the node would occur 
only along k x =0. 



parameter must be substituted into the Eq.(I16>. For simplicity we consider first the solution of Eq. dl6l in 
the Ru02-planes and then discuss what is happening along the c-axis. 

Solving Eq. j!6l > in the first BZ down to 5 K we have found that p-wave symmetry yields the largest 
energy gain for the xy— band, while for the xz— and yz— bands the situation is more complicated. We 
find that the gap equations for the xy- and yz— , rrz-bands can be separated due to their weak interaction. 
The result is that for the xy-band the p-wave is the most stable solution, while for the xz— and yz— bands 
/ 2 ,2_ y 2-wave symmetry yields the largest eigenvalue due to stronger nesting. 

How to understand the competition between p- and /-wave pairing?. For doing this we display in Fig. [5] 
characterizes the solutions of Eq. dl6l . Note, a good approximation is to linearize Eq. J 16b with respect to 
Aj, i.e. E\j — > e l k ,. We use also tanh(ej c ,/2fcsT) = 1. Thus, the main contribution to the pairing comes 
from the states at the Fermi level. We present our results in terms of the Fermi surface of the RUO2 plane. 
The wave vectors Qi and refer to the peaks in x(q, w). The areas with A^ 2 _ 2 > and Af^ 2 _^ 2 < 
are denoted by (+) and (-), respectively. Note, that the minus sign in Eq. H6i is absent for triplet pairing. 
The summation over k' in the first BZ is dominated by the contributions due to Q; for the a- and /3-bands 
and the one due to for the 7-band. 

As can be seen from Fig. |5]in the case of / 2 .2_ J) 2-wave symmetry for the xy— band the wave vector 
bridges the same number of portions of the FS with opposite and equal sign. Therefore, for the xy-band 
no solution with / 2 .2_ J ,2-wave symmetry is possible. On the other hand, for the xz— and yz-band the wave 
vector Qj bridges portions of the FS with equal signs of the f x a _ y 2 -superconducting order parameter. The 
eigenvalue of this order parameter is also enhanced due the interband nesting for xz and yz-bands. Thus, 
superconductivity in the xz- and yz-band in the RuC>2 plane is indeed possible with f x 2_ y 2-wa\e order 
parameter. Therefore, solving Eq. H61 for the three band picture we found a competition between triplet 
p-wave and f x 2_ y 2-wave superconductivity in the RuC>2-plane. 

We can immediately rule out singlet Cooper-pairing. In particular, assuming d x 2_ y 2 -symmetry for 
Sr2RuC<4 the Eq. H6i yields an eigenvalue which is lower than in the case of triplet pairing. As can be 
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seen using Fig. |9]this is plausible. Note, we get for d 2 ,2_j,2-wave symmetry a change of sign of the order 
parameter upon crossing the diagonals of the BZ. According to Eq. (I16> wave vectors around connecting 
areas (+) and (-) contribute constructively to the pairing. Contributions due to and the background 
connecting areas with the same sign subtract from the pairing (see Fig. [9] with nodes at the diagonals). 
Therefore, we find that the four contributions due to q s ; in the xy — band do not allow d x 2_ y 2-wa\e 
symmetry in the xy— band. Despite the pair-building contribution due to Q; one finds that the eigenvalue 
of the d x 2_ y 2-wave symmetry in the xy-band is smaller than for the f x 2_ y 2-wave symmetry. This is due 
to the large contribution from Q; to the cross-terms for the triplet pairing which are absent for singlet 
pairing. For <i xy -symmetry where the nodes are along (tt,0) and (0,7r) directions we argue similarly. Thus, 
we exclude this symmetry. 

The / 2 .2_j / 2-wave symmetry solution has vertical lines of nodes along the fcj-axis. In principle, this 
has to be observable in the strong four-fold anisotropy of the thermal conductivity n(6, H) in applied 
magnetic field along afc-plane. One has to admit, however, the recent measurements have shown only a 
weak anisotropy 1 27 1 that suggests the line nodes in Sr2Ru04 are not vertical. 

A phenomenological mechanism to produce horizontal line nodes has been proposed by Zhitomirskiy 
and Rice 113 01 . As we have already mentioned, the magnetic fluctuations are different with respect to the 
various bands. The ferromagnetic fluctuations are originating from the xy-band while the incommensurate 
antiferromagnetic fluctuations are due to the nesting of xz and yz-bands. In this case it was proposed that 
in the active a;y-band the superconducting gap has no nodes and can be described by the nodeless p-wave 
gap. Then, due to interband scattering of the Cooper-pairs the superconductivity is induced in the xz and 
yz-bands with the symmetry 



Note, this gap has the horizontal line of nodes at k z — ±tt/c. Then even a combination of the nodeless 
p-wave symmetry and that with nodes will still be present, only the position of the nodes will be slightly 
shifted along c-axis. This behavior correctly describes the temperature dependence of the specific heat 
below T c 1 30 1 . A similar model was proposed by Annett et al. 1311 . In this model there are two distinct 
attractive pairing interactions. One, acting between nearest neighbor in-plane xy orbitals, leads to a node- 
less chiral p-wave state on the 7-band. The second acts between nearest neighbors between planes, and is 
assumed to act primarily on the d xz and d yz orbitals, which have lobes perpendicular to the planes. This 
second interaction produces a gap with line nodes on the a and (3 Fermi surface sheets. This gap model 
was shown to be in quite good agreement with experiments on specific heat [25 1, penetration depth |28| 
and thermal conductivity [32 1. Similar to the Zhitomirskiy Rice approach [30 1, this model is also 'orbital 
dependent superconductivity', but in this case there are no 'active' and 'passive' bands, since all bands 
have a pairing interaction and the gap is of similar magnitude on all three Fermi surface sheets. 

Alternatively, interlayer coupling models have been proposed 1 33 34l l35l . These also give rise to gap 
functions of the form of Eq. d23i . either on one or all of the three Fermi surface sheets. The physical mech- 
anisms responsible for the interlayer coupling and hence the gap node are assumed to be either interlayer 
Coulomb interactions, or dipole-dipole interactions. 

The question remains whether spin fluctuations can also explain the formation of the nodes along the 
c-axis in the 'passive' xz and yz-bands. In order to explain this we have first to take into account the 
magnetic anisotropy observed by NMR in the normal state |4|. In the following we will show that spin- 
orbit coupling is an important interaction in the physics of Sr2Ru04 and we shall argue that it leads to the 
formation of the horizontal line nodes in the superconducting state. 




(23) 
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2.4 Spin-orbit coupling effects in the normal and superconducting states 



So far we investigated the three-band Hubbard Hamiltonian. This model neglects spin-orbit coupling and 
would fail to explain the magnetic anisotropy at low temperatures observed in NMR experiments |4|. Re- 
cently it has been proposed that the spin-orbit coupling may play an important role in the superconducting 
state of S^RuOit. In particular, it was shown by Ng and Sigrist [36| that the spin-orbit coupling lowers 
the free-energy of the chiral superconducting order parameter. Another indication of the importance of a 
spin-orbit coupling is the recent observation of the large spin-orbit coupling in the related insulating com- 
pound Ca2RuC>4 1371 . Therefore, we extend the theory by adding to the Hubbard Hamiltonian the spin 
orbit coupling: 



Here, the effective angular momentum 1; operates on the three ^g-orbitals on the site i and are the 
electron spins. As in an earlier approach of Ref. |36| we restrict ourselves to the three orbitals, ignoring 
e2 9 -orbitals and choose the coupling constant A such that the t2 9 -states behave like an Z = 1 angular 
momentum representation. Moreover, it is known that the quasi-two-dimensional xy-band is separated 
from the quasi-one-dimensional xz- and yz-bands. Then, one expects that the effect of spin-orbit coupling 
is relatively small for the xy-band and can be neglected in the first approximation 1381 . Therefore, we 
consider the effect of the spin-orbit coupling on xz- and j/z-bands only. Then, the kinetic part of the 
Hamiltonian H t + H so can be diagonalized and the new energy dispersions are obtained |38|. Note, we 
use A = 80meV in accordance with earlier estimations 1361 1371 . 

Most importantly, the spin-orbit coupling does not break the time-reversal symmetry and therefore 
the Kramers degeneracy between the spin up and down is not removed. The resultant Fermi surface 
consists of three sheets as observed in experiment 1 15 1. Then spin-orbit coupling together with the Hubbard 
Hamiltonian leads to a new quasiparticle which we label with pseudo-spin and pseudo-orbital indices. 
Note, that despite the spin-orbit coupling causing the spin and orbit quantum numbers not to be good ones 
we can still identify the Cooper-pairing to be triplet or singlet. This refers now then to the pseudo-spin 
quantum numbers. At the same time, the magnetic behavior of Sr2Ru04 becomes very anisotropic due to 
the fact that both one-particle Green's functions and Lande g-factors will be different if the magnetic field 
is applied along the c-axis or in the a6-RuC»2 plane. The anisotropy arises mainly from the calculations of 
the Lande g-factors and, in particular, their orbital parts. The factors g z — l z + 2s z and g + = Z + + 2s + are 
calculated using the new quasiparticle states. The latter consist, for example, of xz- and j/z-orbitals which 
in turn are the combinations of the initial orbital states |2, +1 > and |2, — 1 > mixed due to the crystal 
field. Then, the matrix elements < |Z + 1 >, < |Z_ | > are zero for the xz- and j/z-orbitals while the < \l z \ > 
matrix element is not. Therefore, the longitudinal components of the spin susceptibility of the xz- and 
2/2: -band get enhanced in comparison to the transverse one. One of the interesting questions that we will 
analyze later is the effect of spin-orbit coupling on the antiferromagnetic and ferromagnetic fluctuations. 
This provides insights into a microscopic explanation of the pairing mechanism and allows us to calculate 
the spatial structure of the superconducting order parameter. 

For the calculation of the transverse, \i > an d longitudinal, xf z > components of the spin susceptibility 
of each band Z we use the diagrammatic representation shown in Fig. ^3 Note that because the Kramers 

degeneracy is not removed by spin-orbit coupling G l k , = G l k zz . Thus, as we have already mentioned 

before, the anisotropy arises mainly from the calculations of the Lande's g-factors and in particular their 
orbital parts. The factors g z = l z + 2s z and g + = l + + 2s + are calculated using the new quasiparticle 
states. The latter consist, for example, of xz- and yz-orbitals which in turn are the combinations of the 
initial orbital states |2, +1 > and 1 2 , — 1 > mixed due to the crystal field. Then, the same matrix elements 
as above are zero for the xz- and j/z-orbitals while < \l z \ > matrix element is again not. Therefore, the 
longitudinal components of the spin susceptibility of the xz- and yz-band get enhanced in comparison to 
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Fig. 10 Diagrammatic representation of (a) the longitudinal and (b) the transverse magnetic susceptibility. The full 
lines represent the electron Green's function with pseudo-spin a and pseudo-orbital I quantum numbers. The Lande's 
g-factors are denoted by g+ — 1+ + 2s+ (g_ = I- + 2s_) and g z = l z + 2s z . 



the transverse one. We obtain, for example, for the \x z) -states for the transverse susceptibility 



4 x - 2 /( e kxz) /( £ k+qa;z) 

XlxMi") = -^2^(W2k +q -«2kW2k + q) , — ~ — — + 

JV k kxz t k+qxz W lXJ 



(25) 



and for the longitudinal susceptibility 

xo% (q, = xL(q. w ) + xL (q, w) 

^ [«2kU2k+q + «2kf 2k+q + V2(W2kl>2k+q + ^k«2k+q" 



2 

TV 



2 



k 



e te« ~ 4+q^ + W + ZO+ 



(26) 



Here, f(x) is again the Fermi function and u| and u£ are the corresponding coherence factors [38 1 For all 
other orbitals the calculations are similar and straightforward. 

Then, one gets within RPA the following expressions for the transverse susceptibility 

+- ( \ ^r(q^) fryn 
XrpaMM = i-u x +TM ' (27) 

and for the longitudinal susceptibility 

where Xol — Xoi + XoT- These susceptibilities are used in the corresponding pairing interaction for 
triplet pairing. In order to compare our results with NMR and INS experiments we take 

X Z t Z ot = Y,XRPA,i ( 29 ) 



and 



Using the random phase approximation(RPA) for each particular band we calculate the longitudinal and 
transverse components of the total susceptibility in the RuC>2 -plane. Its real parts at u> — are shown in 
Fig. \^\ As a result of the spin-orbit coupling the magnetic response becomes very anisotropic along the 
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Fig. 11 Calculated real part of the longitudinal and transverse components of the total RPA spin susceptibility 
Xrpa* — 12i xt,RPA re f ers 10 the band indices) in the normal state of S^RuCU along the route (0, 0) — > (it, 0) — ► 
(it, 7r) — > (0, 0) in the first BZ. Due to spin-orbit interaction the transverse component (x ) doe not contain incom- 
mensurate antiferromagnetic fluctuations (IAF) from xz and yz-orbitals, while in the longitudinal component (x" z ) 
they are enhanced. 



whole Brillouin Zone. Since the spin and orbital degrees of freedom are mixed now, the orbital anisotropy 
will be reflected in the magnetic susceptibility. As one sees in Fig. ^2 tne longitudinal component has a 
pronounced peak at Q;, while the transverse one does not show these features at all and is almost isotropic 
(it reproduces mostly the response of the rry-band). In order to understand why the longitudinal suscepti- 
bility shows mostly the nesting features of the xz- and yz-bands one has to remember that due to spin-orbit 
coupling the orbital component of the magnetic susceptibility cannot be neglected. Therefore, the matrix 
elements such as < i\l z \j > and < i\l + r_<\\j > have to be taken into account. At the same time the xz- 
and yz-bands consist of the |2, +1 > and |2, — 1 > orbital states. One sees that while the longitudinal com- 
ponent gets an extra term due to l z , the transverse component does not, since the transitions between states 
1 2, +1 > and |2, — 1 > are not allowed. Thus, the contribution of the nesting of the xz and yz orbitals to 
the longitudinal component of the susceptibility is larger than to the transverse one and incommensurate 
antiferromagnetic fluctuations are almost absent. 

In order to demonstrate the temperature dependence of the magnetic anisotropy induced by the spin- 
orbit coupling we display in Fig. ^]the temperature dependence of the quantity 

Ir nxRPA((l,uJ s f) 



for both components. At room temperature both longitudinal and transverse susceptibilities are almost 
identical, since thermal effects wash out the influence of the spin-orbit interaction. With decreasing tem- 
perature the magnetic anisotropy increases and at low temperatures we find the important result that the 
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Fig. 12 Temperature dependence of the magnetic anisotropy reflected by the imaginary part of the spin susceptibility 
divided by u) 3 f and summed over q. Note, zz and +— refer to the out-of-plane (solid curve) and in-plane (dashed 
curve) components of the RPA spin susceptibility. 



out-of-plane component \ zz is about two times larger than the in-plane one: 

X 2Z >X + "/2 • 

We also note that our results are in accordance with earlier estimations made by Ng and Sigrist |39|. 
However, there is one important difference. In their work it was found that the IAF are slightly enhanced 
in the longitudinal components of the xz- and yz-bands in comparison to the transverse one. In our case 
we have found that the longitudinal component of the magnetic susceptibility is strongly enhanced due to 
orbital contributions. Moreover, we have shown by taking into account the correlation effects within RPA 
that the IAF are further enhanced in the z-direction. 

Finally, in order to compare our results with experimental data we calculate the nuclear spin-lattice 
relaxation rate for the 17 ion in the Ru02-plane for different external magnetic field orientation 
(i = a, b, and c): 

= ^ E |^^f£l , (31) 

where is the q-dependent hyperfine-coupling constant perpendicular to the i-direction. 

In Fig. [0]we show the calculated temperature dependence of the spin-lattice relaxation for an external 
magnetic field within and perpendicular to the RuC>2 -plane together with experimental data. At T = 250 K 
the spin-lattice relaxation rate is almost isotropic. Due to the anisotropy in the spin susceptibilities arising 
from spin-orbit coupling the relaxation rates become increasingly different for decreasing temperature. 
The largest anisotropy occurs close to the superconducting transition temperature in good agreement for 
experimental data 1 8 1 . 
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Fig. 13 Calculated normal state temperature dependence of the nuclear spin-lattice relaxation rate T-j -1 of ir O in the 
RUO2 -plane for an external magnetic field applied along c-axis (dashed curve) and along the a6-plane (solid curve) 
using Eq. (31). Triangles are experimental results taken from Ref. 1 8 1. 



Note, most recently similar magnetic anisotropy has been observed by means of inelastic neutron scat- 
tering |40|. However, the absolute magnitude of the anisotropy seems to be smaller than in the case of the 
NMR experiments 1 8 1. The reason for this might be that in the NMR experiments the magnetic suscepti- 
bility is summed over the whole Brillouin Zone while in the INS only its part around Q = (2ir/3, 2tt/3) 
is probed. 

What is the influence of spin-orbit coupling for determining the symmetry of the superconducting gap 
in Sr2RuC>4? First, the spin-orbit coupling affects the spin dynamics and as we have shown induces the 
anisotropy in the spin subspace. In particular the two-dimensional IAF at Qi = (2ir/3, 2tt/3) have a po- 
larization along the z-direction. This simply means that the antiferromagnetic moments associated with 
these fluctuations align along the z-direction. At the same time the ferromagnetic fluctuations are in the 
a6-plane. This is illustrated in Fig. PH] Since the interaction in the RuC>2-plane (k z =Q) is mainly ferro- 
magnetic nodeless p-wave pairing is indeed possible and the superconducting order parameter in all three 
bands will not have nodes. This differs from previous results in which only the xy-band was nodeless. 
Therefore, we safely conclude that triplet p-wave pairing without line-nodes will be realized in the RuC>2- 
plane. On the other hand the state between the Ru02-plane will be determined mainly by the IAF that are 
polarized along the c-direction. This implies that the magnetic interaction between the planes has to be 
rather antiferromagnetic than ferromagnetic. This antiferromagnetic interaction strongly suggests a line of 
nodes for A(k) between the neighboring RuC>2-planes. This is in agreement with the arguments presented 
in Ref. [ 30 1 . Then, due to the magnetic anisotropy induced by spin-orbit coupling a nodeless p— wave 
pairing is possible in the Ru02-plane as experimentally observed 1271 . while a node would lie between the 
Ru02-planes. 

Note, the spin-orbit coupling determines also the orientation of the d-vector of the Cooper-pairs. In a 
triplet superconductor, the pairing state can be represented by the three-dimensional vector d(k), whose 
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Fig. 14 Schematic representation of the possible node formation in the order parameter between the R11O2 -planes (in 
real coordinate representation) as resulting from the magnetic anisotropy. Here, the amplitude of the order parameter 
along the z-direction has been drawn in cylindrical coordinates between RUO2 -planes. This horizontal line node in 
the superconducting gap resulting from antiferromagnetic spin fluctuations along z-direction seems to be in agreement 
with thermal conductivity measurements below T c 1271 . 



magnitude and direction may vary over the Fermi surface in k-space. Like in the case of superfluid 
3 He, several pairing states related to the different orientation of the total spin moment of the Cooper- 
pair d z , d y , d z may have the same condensation energy under weak-coupling conditions. This degeneracy 
is related to the spin rotation symmetry being present in 3 He (which leads to presence of the so-called 
B-phase), but is lifted by the spin-orbit coupling in Sr2Ru04 |36|. Then, the analogy of the A-phase-like 
can be realized in Sr2Ru04 (dCk)^ = (0,0,d z (k x ± ik y ))) while the two-dimensional analog of the _B-phase 
will be suppressed (d(k)s = \d\ (k x , k y , 0)). The different orientation of the spin and orbital angular mo- 
menta of the Cooper-pair is illustrated in Fig. 1151 . The important consequence of this order parameter is 
the uncompensated orbital angular momentum of the Cooper-pair (so-called 'chiral' state). 

In analogy to 3 He, interesting effects can be observed by applying external magnetic field or external 
pressure along a6-plane to Sr2Ru04. For example, in 3 He the transition from the B to the A phase is 
observed under pressure while in an external magnetic field the A transition splits into A\, A2 phases. In 
Sr2RuC<4 the situation is somewhat opposite. Superconductivity in Sr2Ru04 is already in the A-like phase 
and therefore, it would be interesting to see whether by applying pressure or magnetic field one could cause 
a transition into a B-like phase. Note that recently a second phase has been observed in high magnetic field 
applied along a-direction (ff||(100)) with T c ~ 0.5 T c0 ED- The origin of this feature is not clear at the 
moment. 



22 



I, Eremin et al.: Unconventional superconductivity and magnetism in Sf2Ru04 




Fig. 15 Illustration of the possible superconducting states in Sr2Ru04 taken from Ref . 1421 . (a) A (k) = d z (sin k x + 
isinfcj,) and (b) A(k) = d x smk x + d y smk y have different spin and total angular momentum. The A(k) = 
d z (sin k x + i sin k y ) superconducting state has an angular momentum along c-axis (thick arrow) and the spin is lying 
in the Ru02-plane. 

2.5 Comparison of superconductivity in cuprates and Si"2Ru04 

Despite of the similar crystal structures the differences in the electronic states close to the Fermi level 
leads to the different electronic structure of cuprates and ruthenates. As a result the spin fluctuations are 
very different. In both systems the spin fluctuations are strong but while in cuprates the antiferromagnetic 
fluctuations at Q=(7r, tt) in Sr2Ru04 there is a competition between ferromagnetic and incommensurate 
antiferromagnetic fluctuations at Qiaf = (2tt/3, 2tt/3). As we have shown earlier the superconductivity 
in the Sr2Ru04 seems to be induced by mainly the ferromagnetic fluctuations and the role of IAF is in the 
formation of the line nodes in the passive xz and yz-bands between the RuC»2-planes. 

The comparison of the superconductivity in ruthenates and cuprates was done in Ref. 1 43 1 on the basis of 
the t— J — I model. For all cuprate compounds the degeneracy between the copper 3d-orbitals is removed 
by the lattice structure. After some straightforward calculations it can be shown that the hybridized copper 
and oxygen orbitals separate. As shown in Fig. I16f a). the state with highest (lowest) energy has mainly 
d a; 2_j / 2-wave character in the undoped cuprates and the d 9 copper has one missing electron {i.e. a hole) 
and gives the Cu-ion a spin i. Thus, in the absence of doping, the cuprate material is well described by a 
model of mostly localized spin-i states. The other orbitals are occupied and therefore can be neglected. 
However, the Coulomb repulsion between holes in the same orbital is strong and must be included. Due to 
the large Coulomb repulsion the d x 2_ y 2 -orbitals split into two so-called lower and upper Hubbard bands 
(LHB and UHB), respectively as shown in Fig. I16t b). Then, at half-filling the system becomes an insulator. 
Furthermore, due to a very large Ud the splitting of d^.^-band is so large that the oxygen p-band lies 
between the UHB and the LHB. The charge transfer gap in cuprates 6 (S — e p — e^) is smaller than Ud, 
and so underdoped cuprates are charge-transfer insulators. When the cuprates are doped the t — J model 
appears as an effective low-energy Hamiltonian describing the moving of a hole on the antiferromagnetic 
background |44|. 

In contrast to the cuprates, in ruthenates the Fermi level lies in the center of the d xy ,d xz , and d yz 
manifold of states. In the ruthenates considering the active xy-band one could also consider the moving of 
the electrons on the ferromagnetic background. The instability of the system can be analyzed with respect 
to the triplet p-wave and singlet <i-wave Cooper-pairing. In a weak-coupling BCS theory it was found 1 43 1 
that for p-wave pairing the van Hove singularity does not contribute to the T c since in the effective density 
of states the corresponding electron velocity at the Fermi level is cancelled by the angular dependence of 
the p-wave superconducting gap ~ (sin k x ±i sin k y ). At the same time this cancellation does not occur 
for dj;2_ y 2-wave Cooper pairing and thus the van Hove singularity enhances T c . 

We can also contrast the spin fluctuation pairing mechanism in the cuprates and ruthenates. In the 
strong-coupling Eliashberg-like theory of superconductivity in the cuprates 1 45 1 the enhancement of T c 
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Fig. 16 Schematic electronic structure of the Cu02-planes of cuprates reflecting bonding between a Cu 2+ and two 
2 ~ ions, (a) Initial splitting of the bands taking into account crystal field interaction only. We consider only d- 
electrons (holes) of Cu and p x and p y orbitals of the oxygen. The numbers indicate the occupations of the different 
levels in the undoped compound, (b) Role of the on-site Coulomb repulsion and the splitting of 3d a .2„ H 2 -orbital into 
a lower (LHB) and upper Hubbard band (UHB). Ud is the Coulomb repulsion between electrons at Copper. Due to 
8 — e p — ed < Ud the cuprates are Mott-Hubbard insulators without doping. 

occurs due to pronounced antiferromagnetic spin fluctuations and due to a large spectral weight of the 
pairing interaction at high frequencies. Therefore, by comparison of superconductivity in Sr 2 Ru04 and 
high-T c cuprates one may conclude that the higher T c in cuprates can be explained by the stronger anti- 
ferromagnetic fluctuations and the density of states effects (van Hove singularity). This, however, are only 
quantitative arguments since a complete theory of superconductivity in cuprates is still lacking. 

3 Mott-insulator transition in Ca 2 _ x Sr 1 Ru0 4 

Originally, the investigations of ruthenates were concentrated on Sr 2 Ru04 and SrRuC>3, since it was be- 
lieved that in contrast to cuprates in ruthenates only ferromagnetic fluctuations take place. However, it 
was soon realized that the ruthenates reveal much more interesting behavior. In particular, the richness of 
the electronic, magnetic and structural phase transitions in Ca2- K Sr K Ru04 has received a lot of attention 
recently. 

In particular, in C^-^Sr^RuC^ a transition from triplet superconductivity in Sr2Ru04 to an antiferro- 
magnetic Mott-Hubbard insulator in Ca2Ru04 has been found |46|. For various values of x the different 
phase transitions occur as shown in Fig. ^] For example for x < 0.2 this material is an insulator and 
for the interval 0.2 < x < 0.5 C^-^Sr^RuC^ is a metal with short-range antiferromagnetic order. At 
x ~ 0.5 there is a crossover which is accompanied by the sharp enhancement of the ferromagnetic fluc- 
tuations in the uniform spin susceptibility. For x — > 2 the system becomes superconducting with triplet 
Cooper-pairing. 
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Ca x Sr 

Fig. 17 Phase diagram of Ca2- ;c Sr :E RuC)4. P refers to paramagnetic, CAF to canted antiferromagnetic, M to mag- 
netic, SC to the superconducting phase, -M to the metallic phase, and -I to the insulating phase 1 46 1. 

A neutron diffraction study on Ca2- a ;Sr. I ,Ru04 has been made recently |47|, and results in the structural 
phase diagram of this system, shown in Fig. I18l a~). Here, one can see also the experimental data of the 
structural (Tg), magnetic (Tjy) phase transition and the dashed curve refers to the temperature T p where 
the magnetic susceptibility shows a maximum. As one sees for x < 0.2 the metal-insulator transition is 
accompanied by the structural transition as well where both low and high temperature phases have the 
symmetry Pbca with only difference of the parameter c (S-short c, L-long c). In the metallic phase there 
are two structural phase transitions which are influencing the temperature and concentration dependencies 
of the magnetic susceptibility. In Fig. I18f b) we illustrate the evolution of the structural transitions with 
increasing x, starting from Sr 2 RuC>4. The rotation of the RuO^ octahedra around c— axis leads to a phase 
IAi/acd. For x = 1 we have <f> — 10.8. Below x — 0.5 there are addition rotations around the ab-plane 
axis which increase the angle and suppression of the octahedra along the c-axis. 

An attempt to explain the magnetic phase diagram of Sr2_ x Ca x Ru04 has been made on a basis of 'ab- 
initio' band structure calculations |48|, as shown in Fig. \19\ It was shown that the rotations of the RuOg 
octahedra around the c— axis stabilize the ferromagnetic state, since they are sufficient to reduce the pd—ir 
hybridization between the xy-orbital and the oxygen 2p-states. Then the xy-band becomes narrower and 
the van Hove singularity shifts towards the Fermi level. At the same time the xz- and yz-bands are only 
slightly affected. On the contrary, the rotation around the a6-plane changes all three bands completely and 
that increases the nesting of xz- and yz-bands and may stabilize the antiferromagnetic phase |48|. 

In this picture it is easy to understand why the ferromagnetic fluctuation are strongly enhanced around 
the so-called critical doping x ~ 0.5. Using the LDA parameters for the tight-binding energy dispersion we 
calculated the Lindhard response function of the xy-band for Cao.5Si"i.5Ru04 and Sr2Ru04 as shown in 
Fig. |20| One could clearly see that due to vicinity to the van Hove singularity the ferromagnetic response 
of the xy(7)-band is strongly enhanced, at q = 0. Furthermore, since the a and /3-bands are almost 
unchanged the ferromagnetic response gets much stronger than the antiferromagnetic ones indicating the 
dominance of the ferromagnetic fluctuations around this critical point x ~ 0.5. This is in good agreement 
with experiment |47|. However, this model cannot explain the Mott-Hubbard transition in Ca2Ru04 at 
finite temperatures, since Ca2Ru04 remains an insulator even above Tjy. Therefore, strong electronic 
correlations has to be taken into account, as was proposed recently |49|. Here it was suggested that while 
the xz- and yz-bands are split into the lower and upper Hubbard bands, the xy-band is not split even though 
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Fig. 18 Phase diagram of C^-^Sr^RuCU as measured by means of neutron diffraction |47| including the different 
structural and magnetic phases and the occurrence of the maxima in the magnetic spin susceptibility. In the lower part 
the tilting and rotation distortion of the octahedra are shown together with elongation of the basal planes (Ru— small 
points, 0(1)— larger points). All phases are metallic except S— Pbca. 



it is close to the Mott-Hubbard transition. The splitting between xz(a)- and yz(/3)-bands is increased due 
to interband Coulomb repulsion. Indeed, in a mean-field approximation 

U a pn a np -> U a p <n a > np + U a pn a < np > (32) 

the energy dispersions are renormalized 

= e Q + U a p <np > 

£/3 — > £/3 = £/3 + U al3 <n a > (33) 
and the splitting of the a- and /3-bands is given by 

2 

Ae = lp - l a = ep - e a + U aj3 (< n a > - < np >) - t a p + -U a p. (34) 

Here we use the data for the filling of the a and /3-bands from Ref. 1131 . At U a p w 1 eV the splitting is 
Ae ~ 0.8 eV which is comparable with the band width of these bands. Note in this case one will have still 
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Fig. 19 The calculated magnetic phase diagram of Sr2Ru04 including structural distortions, taken from Ref. 1 4-8 1. 
The solid lines are calculated phase boundaries, while the triangles linked by dashed lines correspond to the experi- 
mental data. 
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Fig. 20 The calculated real part of the Lindhard response function, \i0 (<1, 0) along the symmetry points in the BZ in 
comparison for the 7-band for Cao.sSri.sRuCU and Si"2Ru04. Due to vicinity to the van Hove singularity the response 
of the 7-band is strongly enhanced and ferromagnetic (around q = 0). This is in good agreement with experiment 
l47l. 



three Fermi surfaces (as in LDA calculations) due to the xy-band, upper Hubbard band of /3-band and one 
hole-like Fermi surface from the a-band with filling factor « 0.3. 

Due to the rotation of RuC>6-octahedra around the c-axis the width of the xy (7)-band is reduced and 
the system approaches the critical value U ~ U c . As we mentioned for the rotation of the Ru06-octahedra 
around c-axis, the a and /3-bands are not changing in this doping range (0.5 < x < 2.0). Further decreasing 
x and the rotation of the Ru06-octahedra with respect to the afe-plane reduces the bandwidth of all three 
bands so for x = the 7-band is split into the filled lower and empty upper Hubbard bands. A similar 
situation occurs for the /3-band while the a-band is now just completely filled. This allows us to explain 
the formation of the Mott-Hubbard insulator for Ca2RuC>4 in a way consistent with experimental data of 
X-ray Absorption Spectroscopy |50|. Note that while LDA+DMFT results |52| also predict correctly the 
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Fig. 21 Phase diagram of S^Rui-^Ti^C^ 1 57 1. The superconductivity is completely suppressed for x > 0.02 and it 
becomes magnetically ordered with Qi C ~ (2tt/3, 27r/3) indicating an importance of incommensurate antiferromag- 
netic fluctuations of the xz- and yz-bands for the magnetically-ordered state. 

formation of Mott-Hubbard antiferromagnetic insulator in Ca 2 Ru04 the resulting filling of the bands is 
inconsistent with the observed data |50|. 

In addition, a new phase transition in Ca2RuC>4 has been found recently by applying an external hydro- 
static pressure [51]. Above a very small pressure (p c — 0.5 GPa) Ca2RuC>4 shows metallic properties, and 
most unusually it becomes a ferromagnetic metal with a Curie temperature T c strongly dependent on pres- 
sure. For example, just above the critical pressure T c = 12 K and it reaches its maximum value of T c = 25 
K at 5 GPa. In general, in the Hubbard model one would expect a transition from the antiferromagnetic 
insulator to the paramagnetic metal. However, if one keeps in mind the multi-band nature of the metallic 
state, a ferromagnetic ordering is possible. Of course, a study of the mechanism of magnetic order in the 
multiband metallic state in Ca2Ru04 should be done in the future. 



4 Magnetic versus non-magnetic impurities in Sr 2 Ru0 4 

It is known that unconventional superconductors like heavy-fermion compounds or high-T c cuprates reveal 
peculiar behavior if one add magnetic or non-magnetic impurities. In contrast to conventional s-wave su- 
perconductors the non-magnetic as well as the magnetic impurities act as strong pair breakers and suppress 
the transition temperature T c of unconventional superconductors. The suppression of T c reflects the sen- 
sitivity to translational symmetry breaking and is characteristic of anisotropic Cooper-pairing. In cuprates 
the effects induced by magnetic and non-magnetic impurities are not well understood. For example by 
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Fig. 22 Superconducting transition temperature T c as a function of the impurity concentration x (Ti, Ir) 1571 . The 
superconductivity is completely suppressed for x > 0.0015. The dashed curve shows the fit to the Abrikosov-Gor'kov 
pair-breaking function generalized to the case of magnetic and non-magnetic impurities in unconventional supercon- 
ductors |62|. 

doping the non-magnetic Zn the local magnetic moments in Q1O2 -plane are induced around these impuri- 
ties 1 53 1 which demonstrate Kondo-like behavior. On the contrary, the magnetic Ni impurities also show 
quite puzzling features which so far have not been understood [54|. The complication of the cuprates are 
the highly correlated electrons (holes) of the Q1O2 -plane. In this respect it is of high interest to analyze 
the effect of impurities in Sr2Ru04, since it is a Fermi-liquid system in the normal state. 

Early studies of the impurity effects in Si"2Ru04 revealed some features reflecting the unconventional 
nature of superconductivity. For example, the rapid suppression of T c by native impurities and defects 
occurs 1 3 1 as well as the large enhancement of the residual density of states in the superconducting state 
seen in the specific heat [55 1 and NMR measurements 115 61 . 

Quite recently the first accurate study of the effects induced by the magnetic (Ir 4+ ) and non-magnetic 
(Ti 4+ ) impurity substitution in the RuC>2-plane has been performed |57|. Here, the observed effects are also 
quite peculiar. In particular, the substitution of the non-magnetic impurity Ti 4+ (3d ) in S^RuCU induces 
a local magnetic moment with an effective moment p e f f ~ 0.5/j.b /Ti 1581 . The induced moment has Ising 
anisotropy with an easy axis along the c-direction. The corresponding phase diagram is shown in Fig. 
PHI Furthermore, magnetic ordering with glassy behavior appears for x(Ti)> 0.025 in S^Rui-.jTi^C^ 
while retaining the metallic conduction along the in-plane direction. When x(Ti) is further increased 
to x = 0.09 the elastic neutron scattering measurements detect an incommensurate Bragg peak whose 
wave vector Qi C ~ (2tt/3, 2ir/3) is close to the position of the inelastic neutron scattering peak in pure 
Si"2Ru04 115 91 . Most interestingly, in the vicinity of a magnetic ordering a deviation from a pure Fermi- 
liquid behavior seen in Sr2Ru04 is observed by means of resistivity and transport measurements showing 
linear and logarithmic temperature dependencies, respectively |60|. These results indicate that the two- 
dimensional incommensurate antiferromagnetic spin fluctuations arising from the nesting of xz- and yz- 
bands become a static spin density wave state (SDW) by introducing Ti substitution. On the other hand, 
the system Sr2Rui_ a: Ir a ;04 in which the substitution is magnetic Ir 4+ (5c? 5 in the low spin configuration) 
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shows a weak ferromagnetism for x(Ir) > 0.3 161 1 . Thus, substitution of magnetic and non-magnetic 
impurities in Si"2Ru04 leads to a different ground states. 

Despite of these differences both magnetic and non-magnetic impurities act similarly suppressing su- 
perconductivity in Si"2Ru04 for x > 0.0015 |57| as shown in Fig. |22] From this behavior it is clear that 
magnetic and non-magnetic impurities act mainly as a potential scatters for superconductivity in Si"2Ru04 
and that the magnetic scattering does not play any role. Although this behavior could be explained in other 
ways, this observation is also consistent with the existence of a spin triplet state. Magnetic impurities break 
singlet Cooper-pairs mainly because of exchange splitting of the single particle state. The equal spin paired 
state would not be subject to such an effect. 

More detailed theoretical investigation is still needed in order to fully understand the role played by the 
impurities in Sr2Ru04 and illuminate the role played by magnetic fluctuations in the ruthenates. 

5 Conclusions 

In this review we have presented some aspects of superconductivity in Sr2RuC»4 which have not been 
reviewed previously in a unified way so far. Most importantly, we analyzed the role of electronic correla- 
tions and the competition between ferromagnetic and incommensurate antiferromagnetic spin fluctuations 
for the formation of unconventional triplet superconductivity in Si"2RuC»4. We have found that the magnetic 
response is strongly anisotropic and is due to spin-orbit coupling. As we have shown it plays a crucial role 
in the formation of triplet p-wave superconductivity with line nodes between the RuC<2-layers. We fur- 
ther argue that spin-fluctuation-mediated Cooper-pairing scenario in a frame of the three-band Hubbard 
Hamiltonian can explain well the symmetry of the superconducting order parameter in Sr2Ru04. Fur- 
thermore, we point out that strong electronic correlations are responsible for the rich phase diagram of 
Ca2- K Sr 2: Ru04 compound upon changing the Ca content. 

We have also discussed the open issues in Sr2Ru04, like the influence of magnetic and non-magnetic 
impurities on the superconducting and normal state of Sr2Ru04. It is clear that the physics of triplet 
superconductivity in Si"2RuC»4 is still far from being completed and remains to be analyzed more in detail. 

We are grateful to M. Sigrist, M. Braden, and K. H. Bennemann for stimulating discussions. The present 
work is supported by the INTAS Grant No. 01-654. 
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